Aging of the Nonlinear Optical Susceptibility of colloidal solutions 
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Using Z-scan and dynamic light scattering measurements we investigate the nonlinear optics 
response of a colloidal solution undergoing dynamics slowing down with age. We study the high 
optical nonlinearity of an organic dye (Rhodamine B) dispersed in a water-clay (Laponite) solution, 
at different clay concentrations (2.0 wt% - 2.6 wt%), experiencing the gelation process. We determine 
the clay platelets self diffusion coefficient and, by its comparison with the structural relaxation 
time, we conclude that the gelation process proceeds through the structuring of interconnecting clay 
platelets network rather than through clusters growth and aggregation. 



Introduction. Complex systems, including soft- 
materials, are characterized by an exponentially large 
number of metastable states [lj routinely visited while 
the system is evolving towards thermodynamic equi- 
librium if prepared far from equilibrium. This pro- 
cess can extend on timescales much longer than the 
atomic/molecular motions ones, due to the large number 
of metastable states . As a result, the response features, 
like the relaxation time, depend on the specific instant 
{the waiting time t w ) and are determined by the current 
meta-stable state. This mechanism is known as "aging" . 
During aging, the fluctuation-dissipation theory (FDT) 
does not strictly hold, but a generalization has been re- 
cently proposed and numerically tested 0. Accord- 
ing to this generalization the temperature, appearing in 
the FDT as proportionality factor between Response and 
Correlation function, is substituted by and effective tem- 
perature, higher that the actual one. Consequently in 
aging systems, which are subject to fluctuations larger 
than the corresponding equilibrium one, also nonlinear 
effects are expected to be strongly sensible to the pres- 
ence of dynamics slowing down with age. 

In this respect, nonlinear optics experiments in com- 
plex systems can be considered as a promising tool to 
study these out of equilibrium dynamics. However, if on 
one hand the temporal behaviour of the nonlinear op- 
tical coefficients was considered since the beginning of 
nonlinear optics Q and various authors recently reported 
detailed investigations and new models for the time dy- 
namics of specific nonlinear optical effects (like modu- 
lational instability or soliton formation, El U ) • At 
present, only few studies have been devoted to investigate 
the time evolution of non-linear optics properties in soft 
matter. Specifically, in Ref. the frequency dependence 
of X2 has been used to determine the influence of solvents 
on the molecular arrangement of colloidal surface groups 
and how this affects the effective interaction potentials 
controlling gelation. In Ref.s 0, 0], Z-scan measure- 
ments have been used to follow the solvent evaporation 
in thin polymer films. 



Very little is known about the microscopic process 
leading to the arrested state; thus, with the aim of study- 
ing the gelation process in colloidal materials, we use here 
the non-linear optical response, as measured by the Z- 
scan technique, to determine the self diffusion coefficient 
of the colloidal particles, a quantity not directly accessi- 
ble by other techniques. Specifically, in this letter we re- 
port on an experimental investigation of a Laponite clay 
powder dispersed in dye Rhodamine-B aqueous solution, 
displaying both high thermal nonlinear optical response 
and gelation process characteristic of the clay suspension. 
Clear experimental evidence and quantitative assessment 
of aging of the nonlinear susceptibility is reported. 

Sample preparation. Rhodamine B (RhB) solution in 
deionized water (pH=7) was prepared at 0.6 mM con- 
centration. Laponite powder, supplied by Laporte Ltd, 
was then dispersed in RhB solution, stirred vigorously 
until the suspensions were cleared and filtered through 
0.22 fim pore size Millipore filters. We consider the start- 
ing aging time (t w = 0) as the time when the suspension 
was filtered. Four different samples at clay concentra- 
tions between 2.0 wt% and 2.6 wt% were prepared and 
on each of them the two different experiments (DLS and 
Z-scan) described below were contemporarily performed. 

Evidences of aging. In order to provide a clear ev- 
idence of the non stationary dynamics in our samples, 
and for their quantitative characterization, we resorted 
to dynamic light scattering (DLS) measurements, which 
were performed using a homemade correlator [l^ in com- 
bination with a standard optical setup based on a Hc-Ne 
(A=632.8 nm) 10 mw laser, a monomode optical fiber 
and a photomultiplier detector. We observed directly the 
normalized intensity I(q, t) correlation function $2(9, t) = 
(I(q, t)I(q, 0))/(I(q, 0)) 2 , with q the modulus of the scat- 
tering wave vector defined as q — (4irn/X) sin (0/2), be- 
ing the scattering angle — 90° in our experiment. DLS 
data were fitted assuming a correlation function made by 
the sum of an exponential function with relaxation time 
Ti and a stretched exponential function with relaxation 
time T2 and stretching coefficient 0, as reported in 0] 
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and commonly used in photocorrelation measurements 
analysis for liquid-like (ergodic) samples. 

Example of DLS data, together with their best fits, 
are shown in Fig. 0\ as symbols and full lines at the 
indicated waiting times. The mean relaxation time r m 

r m =T 2 0- l T (/T 1 ), (1) 

with T(x) the usual Euler gamma function, can be taken 
as a parameter representing the slow dynamics behaviour 
of the correlation function, while the fast dynamics cor- 
relation time t\ turns out to be constant for our clay 
concentrations during the gelation process, as described 
elsewhere 0. In Fig. 0\ the r m vs. t w data in log- 
linear scale highlight the common exponential growth 
behaviour T m (t w ) = tq exp {^i T t w ) (full line trough the 
data) of the mean relaxation time for different Laponite 
concentrations. These results reflect clearly the aging be- 
haviour of the material, that corresponds to t w dependent 
relaxation times. The specific t w dependence of T m (t w ) 
leads to consider a scaling law, that makes all data col- 
lapse on a single master curve. This curve is reported 
in Fig. 03 where the the DLS derived T m (t w ) data are 
shown as a function of n T t w . 

Linear absorption measurements. In order to show 
that aging has no role in the linear optical response of 
our system, we measured the linear absorption spectra 
by UV-visible spectrometer. In Fig. we report repre- 
sentative absorption spectra of RhB 2.4 x 10 -5 M water 
solution at wt %, 1 wt % and 3 wt % clay concentra- 
tions, showing that the linear absorption of RhB changes 
in presence of Laponite and it can be considered indepen- 
dent on colloid concentration and on the waiting time. 
These samples were prepared following the same proce- 
dure described above, using a lower concentrated RhB 
solution to avoid saturation in absorbance spectra. The 
linear absorbance coefficient ao, for the 0.6 mM Rhb so- 
lution was also found to be ao — 6.2 mm -1 and ao — 4.8 
mm -1 for all our RhB-Laponite dispersions, by measur- 
ing the sample transmission at A — 532 nm confirming 
the absence of aging in the linear optical absorption. 

Aging of the nonlinear susceptibility. Using the Z-scan 
technique [15j we measured both nonlinear absorption 
and refraction: the sample was moved along the z— axis 
through the focal point of the input Gaussian beam, and 
the transmitted power was measured as a function of z in 
the far field using a photodiode behind a small calibrated 
pinhole. In our experiments we used a CW pumped 
diode laser operating at power P=10 mW and wave- 
lenght A=532 nm. The beam was focused by means of a 
75 mm focal length lens giving a 20 /im beam waist radius 
w and a Iq = 8 x 10 6 W/m 2 beam central intensity at 
the focus (z = 0). A photodetector, at distance 310 mm 
from lens focus, was used to probe the light power behind 
a 2 mm aperture and the sample was scanned across the 
focus with a 5 cm translation stage. In Fig. 03 we show 
as an example the Z-scan curves of RhB- 2.2 wt % clay 
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FIG. 1: (A) DLS homodyne correlation functions and (B) Z- 
scan data of Laponite 2.2 wt % dispersion in RhB 0.6 mM 
water solution for different waiting times t w and correspond- 
ing fits (full lines), see text. In (B) is also reported RhB 0.6 
mM scan (open circles). 
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FIG. 2: Mean relaxation time r m vs. waiting time t w (A) 
and vs. the scaled variable u T t w (B) in log-linear scale at 
various clay concentrations; continous lines in (A) are the 
fitted curves. 



dispersion at increasing t w correspondingly to the DLS 
correlation functions (shown in panel A). The dynamics 
slowing down of the system clearly affects the nonlinear 
susceptibility. In Fig. 03 the 0.6 mM Rhb solution (no 
clay) transmittance curve (open circles) is also reported. 
Notice that, similarly to the linear absorption shown in 
Fig. the nonlinear optical response is influenced by the 
presence of the Laponite platlets, while, more interest- 
ingly, contrarily to the linear absorbance, the nonlinear 
susceptibility does actually vary undergoing aging. 

The transmittance curves shape in Fig. 03 corresponds 
to a defocusing nonlinear system (n 2 < 0) as due to 
thermal lens-like behavior, well known in dye solutions. 
Other effects, as for example the electrostrictive one, can 
be ruled out as it would produce a focusing behaviour 
(n 2 > 0). Additionally the expected electrostricitve 
contribution to n 2 for laponite platelets is expected to 
be at least two or three order of magnitudes less than 
those competing to spherical particles 0], as due to 
the smaller particle volume, i.e. many order of mag- 
nitudes of the defocusing n 2 measured in our samples 
(\n 2 \ = 10~ 10 m 2 W _1 ). To fit the experimental data 
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we used the analytical description of Z-scan method re- 
ported in 0, based on the Fresnel-Huygens principle, 
here enriched by including sample nonlinear absorbancc 
effects on transmittance, not considered since now. Pre- 
vious investigations concerning thermal lensing effects, as 
e.g. in Ref. |f8| have shown that the radial temperature 
profile adiabatically follows the beam intensity close to 
the beam centre, while it slowly decays in the tails, due 
to thermal diffusion. In this case the problem can be 
treated perturbatively, by first determining the temper- 
ature profile from the input beam, and then determining 
the beam diffraction out of the sample, as it is the case 
here. As a result the following equations govern laser 
beam phase $(£) and intensity /(£) variations inside a 
nonlinear sample: 



2.0 



d/(C) 
dC 

d£(C) 

dC 



kn(I(Q), 



(2) 



(3) 



where a(io) = a o + c*2^0j with cxq, the linear absorption 
coefficient, and a 2 the nonlinear coefficient, n(/(£)) — 
?io + n 2 I{C,), with no, the linear refractive index, and 
n 2 the Kerr coefficient. In (0 and C is the coordi- 
nate inside the sample along the z-direction and k the 
wavenumber. Using the expression for the transmittance 
reported in 0, by solving equations and for a 
TEMqo electric incident field intensity distribution and 
by determining the electric far field Eq, transmittance 
results as: 



T = 



\E c (z,<S> (z))\< 
\Ec{zM 2 



£(*)r (£(*),.$„(*)) 



(4) 



In Eq. T(x 1 y) is the lower incomplete gamma func- 
tion, 
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1 + z 2 /z 2 Q 
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being z Q = 2kiVg, D the distance between the observa- 
tion point and the beam waist and A the effective sam- 
ple thickness A = [I — exp (~a(I ) L)]/a(I ), with L the 
actual sample thickness. We used as fitting expres- 
sion to estimate both n 2 and a 2 ; the latter is shown in 
Fig-EK vs - tw for the different clay concentrations, while 
from our analysis n 2 values (not reported) are practically 
unchanged during the aging time. Indeed the nonlinear 
refractive index coefficient n 2 is less sensitive to mat- 
ter structural changes because it mainly originates from 
the water density profile arising from the temperature 
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FIG. 3: Absorption spectra of RhB 2.4 1 0" 5 M water solution 
(0 %) and Laponite dispersions in RhB 2.4 1 -5 M solution 
at diffrent clay concentrations and waiting times. 



profile a quantity that only slightly changes during 
aging. This is not the case for a 2 , as the water absorp- 
tion coefficient is negligible and the measured a origi- 
nates only from the Rh-B/Laponite complex. On the 
other hand, from Fig. it is evident that the nonlin- 
ear absorption (a 2 ) is clearly influenced by clay struc- 
tural evolution. According to the theoretical analysis re- 
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FIG. 4: Nonlinear absorption coefficient Q2 vs. waiting time 
t-w (A) and vs. the scaled variable fj, a t w (B) at various clay 
concentrations; continous lines in (A) are the fitted curves. 



ported below, we fitted the Z-scan curves using the func- 
tion a 2 = a + bexp (fj, a t w ), with a, b and /i Q as free 
parameters. The fitting curves are reported in Fig. 0J\ 
as full lines. In Fig. 0)3 the same data collapsed on a sin- 
gle master curve are reported as a function of the scaled 
time t w \i a . The time scaling parameters, /i T and /i Q , ob- 
tained respectively from DLS and Z-scan data, are shown 
in Table [I] ft can be noticed that /i T and fi a values are 
comparable except for the highest concentration, as most 
likely due to the different durations of each experiment 
(i.e. 5 minutes for photocorrelation and about 15 minutes 
for Z-scan single measurement), which are not negligible 
in the presence of a fast gelation process. 

Theoretical analysis. At the microscopic level, the ob- 
served data can be interpreted by taking into account 
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TABLE I: Time scaling parameters 



Clay wt% 


fj, T 


Ma 


2.0 


(1.9 ±0.2) 10" 5 


(2.3 4=0.2) 10 _b 


2.2 


(3.5 4= 0.1) 10" 5 


(3.3 4=0.2) 10" 5 


2.4 


(10.5 4=0.2) 10~ 5 


(10.9 4=0.5) 10~ 5 


2.6 


(28±1)10 -5 


(17±1)10 -5 



the so called Ludwig-Soret effect, |20j i. e. the Laponite 
platelets concentration (c) gradient induced by the pres- 
ence of a temperature gradient. Specifically, as the value 
of Vc depends -at fixed VT- on mass self diffusion 
coefficient, in turn strongly related to the aging time, we 
expect a waiting time dependence of Vc. This mass dif- 
fusion process strongly affects the Rhodamine molecule 
concentration due to their adsorption [2j| on clay disks. 
Indeed, we estimated about 60 Rh-B molecules for each 
Laponite platelet, much less that the maximum allowed 
adsorption sites (about 500 sites as provided by Laporte 
Ltd. datasheet); consequently all dye particles in solu- 
tion can be considered as attached to the colloids. 

The previous argument can be put on a quantitative 
basis by considering the absorption as due to N(t w ,I a ) 
dye particles varying in time due to matter and thermal 
diffusion processes and yielding an overall absorption co- 
efficient 



a(t w ,I a ) = N(t w ,I a ) a(I a ). 



(7) 



Being a = a a 4- a 2 I , expanding N(t w , I ) at first order 
in I 0l and retaining only the terms linear in /„, equation 
(0 reads as 



a(t w ,I ) — N a + 



N a a 2 + a a 



( dN 
\dT 



In 



(8) 



From Eq. (jHJ, the time dependent nonlinear absorption 
coefficient is then written as: 



a 2 (t w ) = N Q a 2 + a Q (dN/dI ) 



(9) 



which undergoes aging because of the t w dependence of 
(dN/dlo), i- e. because of the variation of particles num- 
ber due to diffusive effects. To determine the waiting 
time dependence of dN/dI in 0, we rely on the station- 
ary solutions of both heat equation and mass flow conti- 
nuity equation reported in |20j, which give the following 
radial profile for the concentration c{r,t w ) (ex N(r,t w )): 



V 2 c(r,0 



2 L h Dt 



■c{\-c)e- 2r2 /< (10) 



being k the thermal conductivity, 2 L Iq exp(—2 r 2 / Wq) 
the Gaussian beam source term, c the average colloids 
concentration, D(t w ) the mass self diffusion coefficient 



and Dt the thermal diffusion coefficient. In writing equa- 
tion (TJJJ we assume that the concentration gradient time 
dependence is enclosed uniquely in mass diffusion coef- 
ficient, while Dt was assume to be unchanged, and all 
other terms are constant. Solution of Eq. (|1(J|) yields, 
besides an integration constant, the following radial con- 
centration profile 



c(r, t w )-- 



D T LI w 2 c(l — c) 

*KD(t w ) 



-2r 2 



2log(r) 



(11) 



where Ei{x) denotes the exponential integral function. 
The maximum difference in concentration within the 
sample in a confined geometry of thickness L given by 
dJ is: 



Ac(t w ) 



Dt LIqw 2 c(1 — c) 



[ 7 + 2log(H)}, (12) 



with 7 the Euler constant; in writing i|12|) we regard the 
logaritmic and the exponential terms as main contribu- 
tions to Ac(t w ) for large (r ~ H, being H the transverse 
sample cell dimension) and for small (r ~ 0) distances 
respectively. 

Finally, from Eqs. (O and ©, 



<X2{t w ) = N Q a 2 -C 



D(t w ) 



(13) 



where the constant C encloses all time independent terms 
in 112|) and it also takes account of the proportionality 
relation c(t w ) oc N(t w ). Expression l|13fl gives the t w 
dependence of a 2l and, noticing that /i T w /i Q , one can 
conclude that D oc t^ 1 . This proportionality suggests 
that the gelation process proceeds through the structur- 
ing of interconnecting clay platelets network rather than 
through clusters growth and aggregation. 

In conclusion, we studied the aging of the non-linear 
optical susceptibility in a water-RhB-clay solution. The 
dye molecules are adsorbed on the clay surface, creating 
a complex medium with a high optical nonlinearity ca- 
pable to perform a liquid-gel transition on a time scale 
/i^ 1 , dependent on the clay concentration. To analyze 
the waiting time dependence of the Z-scan curves, we de- 
rived Eqs. (1 IKit . which generalize the result obtained in 
Ref. 01 to the case where photon absorption is present. 
At variance with the linear optical properties, which are 
not affected by aging, we showed that i) the non-linear 
optical susceptibility actually ages, ii) its t w dependence 
is only function of a concentration dependent time scale 
(la 1 , Hi) this time scale is the same as that for the wait- 
ing time dependence of the relaxation time measured by 
DLS, iv) the waiting time dependence can be microscop- 
ically explained assuming a thermodiffusive behavior of 
the clay platelets (Soret effect), and v) the inverse propor- 
tionality between D and r m suggests that the structural 
arrest of the laponite platelets takes place via the orga- 
nization and the structuring of a spatially homogeneous 
clay network. 
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